Physical Review Letters 



Measurements of Surface Diffusivity and Coarsening During Pulsed Laser Deposition 
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Pulsed Laser Deposition (PLD) of homoepitaxial SrTiOs(OOl) was studied with in-situ x-ray 
specular reflectivity and surface diffuse x-ray scattering. Unlike prior reflecivity-based studies, these 
measurements access both the time- and the length-scales of the evolution of the surface morphogy 
during growth. In particular, we show that this technique allows direct measurements of the dif- 
fusivity for both inter- and intra- layer transport. Our results explicitly limit the possible role of 
island break-up, demonstrate the key roles played by nucleation and coarsening in PLD, and place 
an upper bound on the Ehrlich Schwoebel (ES) barrier for downhill diffusion. 

PACS numbers: 68.55.-a, 68.55.A-, 68.47.Gh, 81.15.Fg, 61.05.cf 
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PLD presents an exceptional challenge for experi- 
mental and theoretical study due to its highly non- 
equilibrium nature, the vast range of time- and length- 
scales involved, and the complex stoichiometry of the 
materials system studied. Consequently, fundamental is- 
sues, such as the roles played by the pulsed nature and 
the kinetic energy of the deposit, remain unresolved [3, 0, 
0, 0, HI • System-specific kinetic properties are also diffi- 
cult to obtain. For example, STM has revealed a rich va- 
riety of phenomena on SrTiOs surfaces [(J Q, but at time 
scales longer than those relevant to growth. In contrast, 
fast studies of PLD have typically employed electron @,@] 
or x-ray[3, 0,0,0, [l(| specular reflectivity. These studies 
have excellent time resolution, but are sensitive only to 
the average atomic-scale surface roughness [III [Hj], and 
therefore provide an incomplete description of surface ki- 
netics. 

In this Letter, we show that in-situ x-ray diffuse scat- 
tering provides critical length scale information absent 
from x-ray reflectivity alone, at time scales sufficient to 
study PLD. The experimental details are given in Ref. 
[Hj]. Figs. QJa-c) show false color images of the inten- 
sity of both the specular (q\\ = 0) and the surface diffuse 
scattering as a function of time and q\ \ , during the depo- 
sition of approximately 11 monolayers (ML) of unit cell 
step height at 3 temperatures. As material is deposited 
on the surface, the specular intensity drops while diffuse 
lobes of scattering appear on both sides of the specular 
rod. These lobes are cuts through "Henzler rings" aris- 
ing from 2D islands on the surface 3, lH , as verified by 
ex-situ atomic force microscopy (See Ref. [HI Fig. 2). 

At low layer coverage, (0 < < 0.4 ML), the radius 
of these rings, go, is inversely proportional to the average 
island separation, Li s i ~ 27r/go [ifil]- As more material 
is deposited, the intensity of the specular rod and the 
diffuse lobes oscillate out of phase with a period of 1 ML. 
Near layer completion (0.7 < < 1 ML), qo is a measure 
of the separation between holes rather than islands. 



A conspicuous feature of Figs. QJa-c) is that increas- 
ing the substrate temperature results in a decrease in 
go, corresponding to a decrease in island density, as ex- 
pected from classical nucleation theory [17j. A second 
feature of the data is that qo decreases with increasing 
layer number. This is a general feature of every data 
set we obtained, and reflects the growth surface's "mem- 
ory" of underlying layers. If a new layer nucleates before 
layer completion, the remaining holes function as adatom 
sinks, suppressing the adatom density, thereby producing 
a smaller nucleation density. 

Figs. QJd-f) show an enlarged view of the 1st ML of 
growth from Figs. [TJa-c). At 1000°C, diffuse scattering 
appears between the first and second pulses. At suffi- 
ciently lower temperatures (< 785° C), diffuse scattering 
is not visible until after the second pulse, indicating ei- 
ther delayed nucleation or intensity below our detection 
limit. 

To extract quantitative information, the x-ray data 
were fit to the sum of three independent components, 

I fit{%) = h g +I 8P ec(q\\) +Idiff(q\\ ±g )- W 

In this equation, Ib g is a constant background, and I spec 
and Idiff take the form: 



2„2i3/2 



(2) 



where £ is the correlation length. The parameters 1$ and 
£ each take on two values, associated with I spec and Idiff- 
Eq. ([2]) with qo = corresponds to the scattering profile 
of a random distribution of islands [Til. fl9j |. Figs. OJg-i) 
show the single frames from Figs. QJd-f) corresponding 
to t = 16.5 s: the frame immediately following the third 
laser pulse. Also shown are the best fit to Eq. (pp) and its 
components. The agreement between the fitting function 
and the data is excellent, with a typical x 2 ~ 1-3. 

Fig. [2^l shows the evolution of qo and £ for the first 
monolayer at 850° C. Immediately following the first 
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FIG. 1: Diffuse x-ray scattering for the PLD of SrTiO 3 (001). (a-c) Depositions of -11 ML at 1000°C, 785°C, 695°C, respectively, 
(d-f) The corresponding first ML. Vertical lines represent the laser pulses (first pulse at 5 s). (g-i) Scattering line shape at 
t = 16.5 s for each temperature. If a (red, solid) consists of Idiff (green, dash-dot), I spe c (black, solid) and hg (black, dashed). 



pulse, a diffuse peak is observed at qo = 0.066 ± 0.005 
A, indicating that some islands have nucleated. This 
value of qo corresponds to an island density of n x = 
(1.0 ±0.1) x 10 12 cm -2 if a triangular lattice is assumed. 
A rising qo immediately following the first pulse would 
signify nucleation of new islands from a supersaturation 
of adatoms. Instead, qo decreases monotonically and con- 
tinuously, indicating a steadily decreasing island density. 
This shows that some of the newly formed islands are dis- 
appearing, and thus that island coarsening 20i [2lj . rather 
than nucleation, drives the evolution of qo during this 
time. We observe similar coarsening for substrate tem- 
peratures as low as 695° C. 

A key parameter in PLD growth is the decay time of 
the adatom supersaturation resulting from the pulse [221]. 
Our diffuse scattering measurements are not directly sen- 
sitive to adatom supersaturation. Specifically, since they 
only extend only to q ma x = 0.2 A -1 (see Figs. QJg-i)), 
they are insensitive to lateral correlations smaller than 
~ 2ixjq rnax ~ 31 A, such as adatoms or very small is- 
lands. However, it is easily shown that, if the cover- 
age and specular intensity are both constant, the to- 
tal diffuse scattering intensity is also constant [30] . We 
therefore write the total in-plane surface scattering as 
hot = Lpec + hsi + Lm, where I sm accounts for features 



not captured by our detector. Ii s i is equal to Idiff from 
Eq. (pQ), integrated over the q z plane: 



I isl = 27rJo(«b/0[l + Vl + teZo)- 2 ], 



(3) 



and is associated with the total diffuse scattering due to 
large islands, i.e. the islands separated by > 2i\ j q max . 
When the specular intensity between pulses is constant, 
such as at low coverage in Fig. 2b, a time-dependent Ii S \ 
corresponds to mass transfer between small features and 
the characteristic large islands that give rise to Ii S \ . 

The specular intensity, I spec , and total diffuse intensity 
Iisi for an 850°C deposition are shown in Fig. [2)3. Apart 
from the jumps in I spec associated with each deposition 
pulse, we observe two, distinct slower changes occurring 
between pulses. The first is a change in I spec that oc- 
curs near monolayer completion and has been studied 
previously [I, 0, 0, 0, H[. The second slow change, which 
manifests in Ii s i and has not previously been reported, 
occurs at low coverage. After the third laser pulse, the 
rise in 1^ lags behind the fast drop in I spe c- As discussed 
above, this delay indicates an increase in the amount of 
material in large islands. Moreover, since I sp ec is con- 
stant during this time, this mass transfer corresponds 
solely to intralayer transport. 
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FIG. 2: (a) The peak position of the diffuse lobes, go (A) and 
the correlation length, at 850° C are shown for the 1 st 
ML. Vertical lines represent laser pulses (1 st pulse at 5 s). (b) 
I spec (A) and Idiff (•) are shown. The characteristic diffusion 
times, nsi and Tspec are determined by fitting J- spec and Ii„i. 



The relaxation kinetics described above can be quan- 
tified by fitting I spec at high and Ii s i at low to a 
simple exponential with characteristic relaxation times 
r spec and r^. However, the physical process or pro- 
cesses giving rise to these time constants cannot be de- 
termined from Fig. [2] alone. For example, the diffusing 
species may come from pre-existing islands; therefore, 
both 

T spec and Ti s i may be determined by either the rate 
of adatom detachment or the rate of surface diffusion. If 
present, an ES barrier for downhill diffusion would also 
contribute to r spec . We are able to resolve this ambiguity 
by examining the relationship between r spec and go ob- 
tained for different layers in a single growth, exploiting 
the fact that qo decreases with increasing layer number. 
If diffusion is indeed the rate-limiting process determin- 
ing r spec and if the average diffusion length, L^, is de- 
termined by qo then the Einstein relation, L D = 4Dr 
applies [23[. We associate each go with an approximate 
diffusion length = L is i/2 = 7r/go, (approximately 
half the distance between hole centers), and plot L D vs. 
T spec obtained from approximately the same coverage, 
e.g. ~ 0.8, 1.8, 2.8, for several different layers in FigO 
A clear linear relationship is observed, so that we may 
associate the slope in Fig. [3] with the diffusivity, D. We 
also assign n 3 i to diffusion-limited transport, since only a 
subset of the processes responsible for r spec are involved. 

Figs. [U(a,b) show Arrhenius plots of D obtained from 
the analysis of r spec and ri 3 i for the first ML. The best- 
fit lines are shown, corresponding to activation energies 
of E a = 1.0 ± 0.1 eV and E a = 0.9 ± 0.2 eV for inter- 
and intra-layer transport, respectively. The difference in 
these energies, 0.1 ±0.22 eV, is a direct measure of the ES 
barrier. Remarkably, these data sets yield not only the 
same slopes (within experimental error) but also the same 




FIG. 3: Length scale for diffusion, L 2 D vs. r sp ec for 1000° C (•), 
850° C (o), 785°, (*), and 695° C (□). The linear relationship 
shows that diffusion is the rate limiting process. 



values of diffusivity throughout the temperature range 
studied, suggesting that the ES barrier is negligible. We 
thus combine the data in Figs. S^a-b) to give the single 
result D = D exp(-E a /k B T), with D = 10- 8±1 cm 2 
s -1 and E a = 0.97 ± 0.07 eV. The determination of both 
Do and E a through diffraction-based measurements alone 
represents a principle result of this work. 

The value of E a reported here is larger than two values, 
0.48dz0.05 eV and 0.6±0.2 eV, previously reported in the 
liter ature [2|, ll0[. In these reports, E a was obtained from 
the temperature dependence of r spec implicitly assuming 
a constant length scale. The effect of this assumption 
on the determination of E a is made explicit by writing 
the temperature dependence of the length scale in Arrhe- 
nius form, Lb = Lq exp(— El/J^sT), and rewriting the 
Einstein relation: 



rspec = (L 2 /4D )exp[(E a - 2E L )/k B T] 



(4) 



Eq. (|4|) shows that the activation energy measured 
from T spec alone underestimates the activation barrier for 
diffusion, E ai by 2E B . We note that our value of E a = 
0.97±0.07eV is very close to that of 1.2±0.1 eV measured 
for diffusion of TiO x "diline" units on a reconstructed 
SrTi0 3 surface 0. 

Our results provide new insight into the possibility 
of energetic mechanisms promoting smooth growth in 
complex oxide PLD. One such proposed mechanism is 
island breakup, in which energetic impinging material 
breaks up existing islands, delaying second-layer nucle- 
ation. Island-breakup has previously been observed in 
simulations of metal/metal epitaxy [3, HI] and was re- 
cently invoked [5] to explain experimental results of PLD 
of Lai_^Sr^Mn03 on SrTiOs. Specifically, Ref.[5] sug- 
gests that island break-up produces an increasing island 
density when 6 < 0.5 ML. Although the system stud- 
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FIG. 4: (a) Arrhenius behavior of the diffusivity at ~ 
0.8ML. (inset) r spec is obtained from the specular relaxation 
at high coverage, during inter- layer transport, (b). Diffusiv- 
ity at at « 0.25 ML. (inset) Ti S \ is determined by fitting the 
time evolution of U s i. 

ied here is not precisely the same as in Ref.[5], Fig. [2^ 
demonstrates that the island density monotonically de- 
creases with 6 from the earliest moments after nucleation. 
Island breakup could also manifest in our measurement 
as a decrease in 1^ as mass is transferred from large is- 
lands to smaller species without changing qo. However, 
we do not observe such a decrease. Thus, island coars- 
ening overwhelms any possible effect of island break-up. 

A second proposed non-thermal smoothing mechanism 
suggested by prior experimental work on complex-oxide 
PLD, is enhanced downhill transport 0, 0, H Hi. The 
experimental basis for this suggestion is the observation, 
based on specular reflectivity, that downhill transport oc- 
curs on two widely separated time scales ^ Our ob- 
servation, that island nucleation occurs quickly, followed 
by coarsening, suggests an alternate origin of these two 
time scales. Specifically, it is possible that the mobile 
species responsible for slow downhill transport consists 
of material that detaches from islands. This material 
need not be chemically identical with the species arriving 
from the plume. Interestingly, we note that the prefac- 
tor reported here, Do = 10 _8±1 cm 2 s _1 , is five orders of 
magnitude lower than typical experimental and theoreti- 
cal value for metal and semiconductor systems [26]. Sim- 
ilar diminished prefactors have previously been associ- 
ated with correlated motion involving multiple atoms [26j. 
Here, it might be associated with stoichiometric mass 



transfer of Sr-containing and Ti-containing species. 

In summary, we have presented time-resolved x-ray re- 
flectivity and diffuse scattering measurements obtained 
during PLD. Our results constitute direct observations of 
island nucleation as little as 200 ms after the pulse, and 
direct evidence of island coarsening occurring between 
laser pulses for temperatures as low as 695° C. Quanti- 
tative analysis of our results allow us to independently 
estimate the inter- and intra-layer diffusivity (prefactor 
and activation barrier) of mobile species between pulses 
and to place an upper bound on the ES barrier. Our 
measurements significantly impact prior estimates of the 
thermal diffusivity involved in SrTi03 growth, and place 
specific constraints on energetic smoothing mechanisms 
that have been proposed to occur during complex oxide 
PLD. 
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EXPERIMENTAL DETAILS. 

The x-ray measurements were performed using a cus- 
tom PLD/x-ray diffraction system installed in the G3 
hutch at the Cornell High Energy Synchrotron Source. 
Auxiliary Figj5] illustrates the experimental geometry. 
The surface shown is an atomic force microscopy (AFM) 
image taken after the deposition performed in Fig. 1(a) 
of the manuscript where 6 = 11.3 ML. As suggested by 
the figure, steps on the sample resulting from the surface 
miscut are always aligned perpendicular to the incident 
beam. The x-ray beam produces scattering in the qu 
direction that is directly related to the island correla- 
tions on the surface. For these experiments, a monochro- 
matic (AE/E=1%) 10.0 keV x-ray beam with 8 x 10 13 
photons/sec/mm 2 was slit down to produce a 1.0 mm x 
0.5 mm beam at the sample. The specular component of 
the reflected signal was attenuated using a tin absorber 
to prevent detector saturation. To optimize the signal- 
to-background ratio of the diffuse scattering, the exper- 
iment was performed near the L = 0.275 r.l.u. position 
on the crystal truncation rod. The in-plane diffuse sur- 
face scattering was monitored using a CCD area detector 
operating as a linear detector in streak mode. The time 
resolution of the experiment is limited by both the read- 
out time of the detector 78 ms) and by the incident 
x-ray flux (100 — 200 ms). For these growth conditions, 
we are able to collect 18 images between laser pulses. 

The depositions were performed by laser ablating a 
single crystal SrTiOs target using a 100 MW/cm 2 KrF 
excimer laser (248 nm). The target is located 6 cm from 
the substrate. The area of the laser spot on the tar- 
get was approximately 3.7 mm 2 with a fluence of 1.9 
J/cm 2 . This configuration deposited ~0.09 ML/pulse at 



FIG. 5: Diffuse Scattering Diagram. A 10 keV x-ray beam 
reflects off the surface with a specular (purple) and diffuse 
scattering (gold) which is detected with a CCD detector. The 
value of q\\ measured is directly related to the island correla- 
tions on the surface. 



a laser repetition rate of 0.2 Hz, with a 2 x 10 -4 Torr 
partial pressure of O2. The substrate temperature was 
measured using an optical pyrometer (A = 4.8 — 5.3/im, 
emissivity=0.8). The substrate preparation procedure 
employed [2 71 l2£ 



FIG. 6: AFM after 11.3 ML Deposition: Post deposition AFM 
of the 1000C growth reveals unit cell high islands on the sur- 
face. The diffuse ring present in the fast fourier transform 
(FFT) shows that the islands are correlated (inset). 



produced a Ti02 terminated surface, 
and AFM confirmed the presence of single unit cell high 
steps separating large atomically flat terraces. 

The diffuse scattering peaks are the direct result of 
single unit cell high islands, which was confirmed by ex- 
situ AFM. An AFM image after the deposition at 1000C, 
shown in Fig. 1(a) of the manuscript, is given in Auxil- 
iary Figj6] confirming the presence of islands on the sur- 
face. The fast Fourier transform (FFT) of the image is 
shown in the inset. The presence of the diffuse Henzler 
ring is a direct result of the correlated islands on the 
surface. Our experimental detected x-ray intensity is a 
cut through these Henzler rings. Coarsening, during the 
post deposition cool-down (« 3 hours), has caused the in- 
crease of the measured length scale of the AFM relative 
to the x-ray data. 
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[30] By integrating Eq. (6) in Ref. [29[ over q\\, it is seen that transform of the vertical height distribution, 

both I spec and Idiff depend only on h(q z ), the Fourier 



